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We investigate the magnetic field dependence of the plasma resonance frequency in vortex crystal state. We 
find that low magnetic field induces a small correction to the plasma frequency proportional to the field. The 
slope of this linear field dependence is directly related to the average distance between the pancake vortices in the 
neighboring layers, the wandering length. This length is determined by both Josephson and magnetic couplings 
between layers. At higher fields the plasma frequency is suppressed collectively and is determined by elastic energy 
of the vortex lattice. Analyzing experimental data, we find that (i) the wandering length becomes comparable 
with the London penetration depth near T c , (ii) at small melting fields (< 20 G) the wandering length does not 
change noticeably at the melting transition demonstrating existence of the line liquid phase in this field range, and 
(iii) the self consistent theory of pancake fluctuations describes very well the field dependence of the Josephson 
plasma resonance frequency up to the melting point. 



1. Introduction 

Josephson coupling characterizes the ability of 
layered superconductors to carry supercurrents 
across the layers. In very anisotropic supercon- 
ductors this coupling is suppressed by magnetic 
field applied along the c-axis. Thermal fluctua- 
tions and uncorrelated pinning lead to misalign- 
ment of pancake vortices induced by the mag- 
netic field || (see Fig. |l|). Misalignment results in 
nonzero phase difference and in the suppression of 
the Josephson interlayer coupling. This suppres- 
sion is quantitatively characterized by the "lo- 
cal coherence factor" C = (cos (r)), where 
Vn,n+i( r ) is the gauge-invariant phase difference 
between layers n and n + 1, (...) means aver- 
age over thermal disorder and pinning. Josephson 
plasma resonance (JPR) measurements in highly 
anisotropic layered superconductors Q probe di- 
rectly the interlayer Josephson coupling and the 
effect of pancake vortices on this coupling, be- 
cause the squared JPR frequency, lj^, in the most 
part of vortex phase diagram, is proportional to 



*This work was supported by the NSF Office of the Sci- 
ence and Technology Center under contract No. DMR-91- 
20000 and by the U. S. DOE, BES-Materials Sciences, un- 
der contract No. W-31- 109-ENG-38. Work in Los Alamos 
is supported by U.S. the DOE. 



the average interlayer Josephson energy 
u? ulC oc J C, 



(1) 



where ojq{T) = c/ y/eo\ c (T) is the zero field 
plasma frequency, A C (T) is the c-component of the 
London penetration depth, eo is dielectric con- 
stant, and Jq is the Josephson critical current. 




Figure 1. Meandering of pancakes along the vor- 
tex line in layered superconductors and grayscale 
plot of cos </> rain+ i near two misaligned pancakes. 

The JPR measurements performed in the liq- 
uid vortex phase at relatively high magnetic 
fields, B > Bj = <5> /Aj, revealed that the 
plasma frequency drops approximately as \j\f~B 
||, where Aj = 7s is the Josephson length, 7 
is the anisotropy ratio and s is the interlayer dis- 
tance. The above dependence is characteristic for 
the pancake liquid weakly correlated along the c 
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axis.y In this phase many pancake vortices con- 
tribute to the suppression of the phase difference 
at a given point. In contrast, in the vortex solid 
pancake vortices form aligned stacks and suppres- 
sion of coupling is caused by weak misalignments 
of the pancake vortices due to the thermal fluctu- 
ations and random pinning. JPR measurements 
in Bi 2 Sr 2 CaCu 2 8 _ 5 (Bi-2212) crystals §,§ have 
shown that the JPR frequency decreases approx- 
imately linearly with field in the vortex solid. In 
the fields above 20 Oe the interlayer phase co- 
herence changes drastically at the transition line, 
implying the decoupling nature of the first-order 
melting transition in agreement with theoretical 
expectations (see, e.g., Ref. Q). On the other 
hand, at smaller field phase coherence does not 
change considerably at the melting point. || 

In this paper we consider the Josephson cou- 
pling and JPR in the vortex lattice. We focus on 
the suppression of coupling due to thermal fluc- 
tuations of pancake vortices near the equilibrium 
crystal positions, and neglect influence of pinning 
potential. For real Bi-2212 crystals this approx- 
imation is justified at sufficiently high tempera- 
tures (> 40 K). This problem has been considered 
in the past in the simple limiting cases. |]J| How- 
ever quantitative calculation suitable for compar- 
ison with existing JPR data ||[| in a wide field 
range is still absent. At small fields, when vor- 
tices act independently, uj 2 decreases linear with 
B. The linear dependence was observed exper- 
imentally in Refs. ]5||| in solid state in Bi-2212 
crystals. In fields below 20 Oe near T c this linear 
dependence extends to the liquid state providing 
evidence for a line structure of the vortices in the 
liquid at low fields. The regime of independent 
vortices has been considered in Ref. ||. In this 
paper we extend our consideration to higher fields 
up to the melting field. 

2. Low fields. Isolated vortex lines 

Consider small magnetic fields B <C Bj,B\ = 
$o/47rA^ b ,. At these fields regions of suppressed 
coupling are localized near the vortex lines (pan- 
cake stacks) and do not overlap (single vortex 
regime) || . The field- induced change in C in this 
regime is given by SC = 1 — C = BI/$q, where 
I = f d 2 r (1 - cos (<£„, n+ i (r))), and ^„, n+ i(r) is 



phase difference induced by fluctuation displace- 
ments u„ in a single line. The same integral deter- 
mines the tilt stiffness due to the Josephson cou- 
pling. We split integration domain in / into two 
region, r < R and r > R, where R is the interme- 
diate scale r w < R < 7s, with = (u^„ +1 ) and 
uv^n+i = u n+ i — u n . At r < R we can neglect 
screening due to the Josephson currents and take 

/ s , y-Uyn+i , y-u yn 

(fin.n+i (r) = arctan arctan — . 

At r > R we can take y>„ in +iin linear approxima- 
tion with respect to u„, 

(p n ,n+l (r, k z ) ^7^k z [u(fc 2 ) X V] K (^~J > 

where u(k z ) = s J2 n exp(~isk z n)u n , k z = 
(2/s) sin(sfc z /2), and Ko(z) is the modified Bessel 
function. At intermediate distances r w <C r <C 
A,/, both expressions give the same simple result 
ip n ,n+i(r) « [r x u„,„ + i]/r 2 . Using above asymp- 
totics of (fn^n+i we obtain 

I= \ J ^ (!-««(***)) KM 2 

(u£ >n+ i(l-cos(sfc z ))) ' 

Weak logarithmic dependence on u nj „ + i leads to 
the nonharmonic tilt energy. || In the following 
we will eliminate this nonharmonicity using the 
self consistent harmonic approximation (SCHA), 
which results in the substitution ln(A/w 2 „ +1 ) — > 
ln(0.24^4/r^). Approximate evaluation of the 
above integral gives a simple practical relation 
connecting the field-induced correction of the 
plasma frequency aj p (B,T) with r w for the case 
T w < Xj < a: 

wl(T)-wl(B,T) _ nBrj Q.SAj 

wg(T) 2$ n r w U 

This relation allows one to extract r^, from the 
plasma resonance measurements. 

We now calculate when wandering of the 
vortex lines is caused by thermal fluctuations. In 
the single vortex regime rj, is determined by the 
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wandering energy consisting of the Josephson and 
magnetic contributions, 



dk z 
2ir 



sijk, + w M 



\u{k z )y 



(3) 



where e±j ~ (eo/7 2 ) hi (l.33j / (r w k z fj is the line 
tension due to the Josephson coupling, wm ~ 
( £ o/A^ b ) ln(1.5A a 6/r w ) is the effective cage po- 
tential, which appears due to nonlocal magnetic 
interactions between pancake vortices in differ- 
ent layers (it describes the magnetic tilt stiff- 
ness at wave vectors k z > l/r w )pO|,^|, and e = 
<f> 2 /(47rA a b) 2 . Assuming Gaussian fluctuations of 
pancake vortices we obtain 



8T 



f 



sw M i + c + \/r+T 



(4) 



where the parameter ((T) w 4A 2 h (T)/A 2 7 de- 
scribes the relative roles of the Josephson and 
magnetic interactions. Substituting this result 
into Eq. (0) we obtain 



^\ 2 ab BT 
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(5) 



This result of the single vortex regime is valid 
in both solid and liquid states for B <C Bj, be- 
cause in this field range wandering of lines at 
short scales does not change much at the melt- 
ing point. Eq. (||) describes fairly well the sup- 
pression of the plasma frequency at small fields 
js) . With the data of Ref. S for underdoped Bi- 
221 2 with T c « 84.5 K Eq. (§) gives unexpectedly 
large wandering length r w « 1/im at 77 K, which 
is comparable with both X a b and Xj at this tem- 
perature. However, we found that this estimate 
is a in good agreement with the theoretical cal- 
culation (||). 

3. High fields, B > Bj, B\ 

As the field increases two competing effects 
start to influence pancake fluctuations and field 
dependence of the average Josephson energy. 
Vortex interactions diminish pancake fluctua- 
tions. On the other hand, collective suppression 
of the Josephson energy decreases tilt stiffness 
and enhances pancake fluctuations. Using general 



relations connecting the phase perturbations with 
the elastic lattice deformations (see, e.g. Ref. §]]), 
we obtain 



SC 



[<Pn,n+l] 



(27rsn v Y 



2 2 
x fd 2 qdk z ^ ^([Qxq] 2 (u^(^Qq) 2 K 



(6) 



><q,r 



q 2 (q + Qr + fc z 2 /7 2 



where n v = B/&o, Q are the reciprocal lattice 
vectors (the cut off q rn ~ 2.2/r w in the summa- 
tion over Q is established by comparison with the 
single vortex regime), and u;(q, k z ) (u t (q, k z )) 
are the longitudinal (transversal) elastic lattice 
displacements. Integration with respect to the 
in-plane wave vector q is limited by the first Bril- 
louin zone, which we approximate by the circle 
q < K , Kq = 4irn v . For thermal Gaussian 
fluctuations the mean squared averages (u 2 ; ) = 
(|u t ;(q, k z )\ 2 ) are determined by the correspond- 
ing components of the elastic matrix $ t ;(q, k z ) 

(|u M (q,fc z )| 2 )=T/$ M (q,M- 
At high fields B > Bj, B\ : and at large k z , k z ^S> 
1/A b, ^/nZ we have 

$t{q, k z ) = C 66 q 2 + $44(9, k z ), 
<S>i(q 7 k z ) = $n(q) + $u(q,kz), 



1 



is the compression 



where $n « |£ (l 

stiffness, Cq§ = A§§n v eo/A is the shear mod- 
ulus, parameter A§§ < 1 describes fluctuation 
suppression of Cqs, which we approximate as 
Am = 1 — 0AB/B m with B m being the melting 
field, and 



*44(g, k z 



n v e k 2 z 



ln- 



O.lla 2 



n v e 



2A 2 b 
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27 2 rl (l-0.53q 2 /K 2 y 
0.13a 2 \ B 2 k 2 



4nX 2 ab fc|+ 7 V' 

is the tilt stiffness, computed within SCHA, with 
a = l/^/n^. The wandering length r w has to be 
determined self-consistently from the equation 

f d 2 qdk z 



= 2 



(2tt, 

x «N(q,fc z )| 



(1 



cos(sk z )) 
+ (|u*(q,fc 2 )| 2 )). 



(7) 
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Expression (g) for SC can be naturally split into 
the collective contribution SC co u, corresponding 
to Q — term in the Q-summation, 

(2 7 rsn 1 ,) 2 /- dVfc z k 2 z q 2 (u 2 ) / q \ 

^C co/; / —73-- ; -2 , (») 

2 ( 27r ) (q 2 + fc2/ 7 2j 

and the local contribution (SC/ oc coming from Q > 
terms. At high fields B » $oAj we obtain 
approximate expression for SCi oc , which resembles 
the single- vortex result (||) 



in itn v r 2 w 0.58a 
oCioc w — - — In . 



(9) 



If we consider suppression of coupling in the cylin- 
drical Bravais cell near the chosen vortex line, 
then the local term determines suppression of 
coupling caused by this vortex line and the col- 
lective term determines suppression of coupling 
caused by all other vortex lines. In general, rel- 
ative role of the collective term in SC grows with 
field. We use above expressions to calculate the 
field dependence of C for comparison with JPR 
data. 

Recently detailed measurements of field depen- 
dence of JPR frequency in the vortex crystal state 
of Bi-2212 have been done by M. Gaifullin et al 
using frequency scan.|| To compare our calcula- 
tions with JPR data we need to know X a b(T) and 
7 = Xc/^ab- A c is extracted directly from JPR 
frequency at B = 0, A C (T) = c/^/eotoo(T) taking 
eo = ll and 7 is chosen as a fitting parameter. 
Fig. H compares the computed dependence C(B) 
with (uj p (B)/ujp(0)) 2 for three values of temper- 
ature. We also show obtained values of A = A a (, 
and 7. We obtain 7, that slightly grows with tem- 
perature (from 460 at 40 K to 510 at 60 K). Such 
enhancement of 7 is expected due to the phase 
fluctuations. 

In conclusion, we have calculated the field de- 
pendence of the JPR frequency in the vortex crys- 
tal. In the single vortex regime at low magnetic 
fields the JPR provides a direct probe for mean- 
dering of individual lines. The wandering length 
extracted from the JPR data is in agreement with 
the theoretical calculations. Our theory of pan- 
cake fluctuations gives a very good description of 
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Figure 2. Comparison of the field dependence of 
the reduced plasma frequency squared (courtesy 
of M. Gaifullin and Y. Matsuda Q) and calcu- 
lated field dependence of the local coherence fac- 
tor C. 
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